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Abstract The 3C-like protease (3CL pro ) of severe 
acute respiratory syndrome associated coronavirus 
(SARS-CoV) is vital for SARS-CoV replication and is 
a promising drug target. Recombinant 3CL pi ° was 
expressed in Pichia pastoris GS115 as a 42 kDa 
protein that displayed a K m of 15 db 2 pM with 
Dabcyl-KTSAVLQSGFRKME-Edans as substrate. 
Purified 3CL pi ° was used for inhibition and kinetic 
assays with seven flavonoid compounds. The IC 50 of 
six flavonoid compounds were 47-381 pM. Quercetin, 
epigallocatechin gallate and gallocatechin gallate 
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(GCG) displayed good inhibition toward 3CL pi ° with 
IC 50 values of 73, 73 and 47 pM, respectively. GCG 
showed a competitive inhibition pattern with K t value 
of 25 zb 1.7 pM. In molecular docking experiments, 
GCG displayed a binding energy of —14 kcal mol -1 to 
the active site of 3CL pi ° and the galloyl moiety at 3-OH 
position was required for 3CL pi ° inhibition activity. 
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Introduction 

In 2002, the first reported outbreak of severe acute 
respiratory syndrome (SARS) occurred in Guangdong 
province, China. It rapidly spread to over 32 countries 
in Asia, North America, and Europe. The mortality 
rate was approximately 10%, according to World 
Health Organization data. The causative agent of 
SARS is a novel human coronavirus (CoV), desig¬ 
nated SARS-CoV. An efficient therapy and a vaccine 
are not currently available (Xu et al. 2005). SARS- 
CoV is an enveloped positive single-stranded RNA 
virus (Rota et al. 2003) that encodes two proteases for 
proteolytic processing: a papain-like cysteine protease 
(PLP2 pi °) and a chymotrypsin-like cysteine protease 
(3C-like protease; 3CL pi °) located in the non-struc¬ 
tural protein regions nsp3 and nsp5, respectively. 
Since 3CL pi ° is essential for the viral life cycle, it is an 
attractive target for the development of antiviral drugs 
directed against SARS-CoV and other CoV infections 
(Grum-Tokars et al. 2008). 

Flavonoids are a large group of naturally occurring 
phenolic compounds ubiquitously distributed in the 
plant kingdom. Over 4,000 varieties of flavonoids have 
been identified (de Groot and Rauen 1998). Flavonoids 
can act as potent antioxidants, and display anti-inflam¬ 
matory, antiallergic, antihemorrhagic, antimutagenic, 
antineoplastic, and hepatoprotective activities (Tapas 
et al. 2008; de Groot and Rauen 1998). In addition, 
flavonoids inhibit the catalytic activities of a great 
variety of enzymes, including hexokinase, phospholi¬ 
pase C, protein kinase C, a-glucosidase, and a-amylase. 
(Tadera et al. 2006). Some flavonoid compounds such as 
quercetin, quercetin derivatives, catechin, epicatechin, 
epicatechin gallate and epigallocatechin gallate inhibit 
SARS-3CF pi ° expressed in Escherichia coli (Chen et al. 
2005, 2006; Yi et al. 2004). However, there has been no 
report on the high inhibition activity of gallocatechin 
gallate (GCG) against SARS-3CF pro , or the structure- 
activity relationship activity among the aforementioned 
flavonoid compounds. 

Herein, we report on the expression of SARS- 
3CF pi ° in Pichia pastoris GS115, its’ in vitro inhibi¬ 
tion by seven flavonoid compounds belonging to four 
groups of flavonoids (flavonol, flavanonol, isoflavone, 
and flavan-3-ol), and the structure-activity relation¬ 
ship between them. The detailed mechanism of GCG 
inhibition was investigated by enzyme kinetic and 
molecular docking studies. 


Materials and methods 

Preparation of recombinant 3CF pro 

The procedures for construction, transformation, and 
screening for the catalytic domain of SARS 3CF pi ° 
were according to the manufacturer’s instructions 
(Invitrogen, Carlsbad, CA, USA). The gene encoding 
SARS 3CF pi ° polyprotein (amino acid residues 
3,241-3,546, GenBank accession no. AY274119) 
(Benson et al. 2011) was optimized by replacing rare 
codons with high-frequency codons, which were 
selected on the basis of the findings for codons usage 
in Pichia pastoris (Supplementary Fig. 1), and were 
synthesized and cloned in the pUC57 vector (pUC57- 
3CF) by a custom gene synthesis service (GenScript, 
Piscataway, NJ, USA). The 3CF pi ° gene was isolated 
from the vector pUC57-3CF by cutting with EcoRl 
and Notl, and was subcloned into the EcoRl/Notl 
restriction sites of the expression vector pPICZaA. 
The novel construct was named pPICZaA-3CF pi ° and 
transformed into E. coli DH5a (Promega, Madison, 
WI, USA) using a standard heat shock method. 
Transformants harbouring pPICZaA-3CF pro were 
selected from FB agar low salt medium [1% (w/v) 
Tryptone, 0.5% (w/v) yeast extract, 0.5% NaCl, pH 
7.5] containing 25 pg Zeocin ml -1 (Invitrogen). 
Plasmid pPICZaA-3CF pro was amplified in E. coli 
DH5a, linearized by Sacl digestion, and transformed 
into P. pastoris GS115 by a modified FiCl method. 
pPICZaA vector was also linearized by Sacl digestion 
and transformed into P. pastoris GS115 as a negative 
control strain. Screening for positive clones was done 
by PCR with two sets of primers (Bioneer, Deajeon, 
Korea): set one contained 3CF pi ° primers [3CF-F 
(5 / -GTGGATTCAGAAAAATGGCC-3 / ) and 3CF-R 
(5 / -CCGCCTGAAAAGTAACTCCT-3 / )] and set two 
contained a-factor and 3AOX1 primers. The PCR 
conditions were 94°C for 5 min, followed by 
25 cycles of 94°C for 1 min, 53°C for 30 s, 72°C for 
1 min, and a final step of 72°C for 5 min. PCR was 
conducted using PCRmix and the PCR product was 
analyzed by agarose gel electrophoresis. 

Recombinant 3CF pi ° was expressed according to 
the manufacturer’s instructions (Invitrogen). Farge- 
scale expression of 3CF pro was carried out in a 10 1 
fermenter with 4 1 BMMY medium. The fermentation 
conditions were 28°C, pH 6.0, 350 rpm, and aeration at 
1 vvm. Methanol was added as described above. 
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Protein expression by pPICZaA transformed into 
P. past oris GS115 was included as a negative control. 
The yeast was separated from the broth by centrifu¬ 
gation at 8,000 xg for 15 min. The pellet was 
discarded and the supernatant was exchanged by 
20 mM Tris/HCl buffer (pH 7.5) using a Millipore 
membrane. The supernatant was used for ammonium 
sulphate fractionation (from 0 to 85%). The obtained 
proteins were dissolved in 20 mM Tris/HCl buffer (pH 

7.5) and dialyzed against 20 mM Tris/HCl buffer (pH 

7.5) . The detection of recombinant 3CL pi ° in the 
culture supernatant was by 12% SDS-PAGE and 
western blot of the electrotransferred proteins accord¬ 
ing to to the manufacturer’s instructions (GE Health¬ 
care, Buckinghamshire, UK). Activity of the 
recombinant enzyme was detected as described below. 

The proteolytic activity of 3CL pro was measured 
using a fluorescence resonance energy transfer 
(FRET)-based assay with a substrate labeled with 
5 - [(2 / -aminoethyl)-amino] naphthelenesulfonic acid 
(Edans) and 4-[{4-(dimehtylamino)phenyl}azo]ben- 
zoic acid (Dabcyl) as the energy transfer pair 
(Bachem, Bubendorf, Switzerland). The Dabcyl- 
KTSAVLQSGFRKME-Edans fluorogenic peptide 
was used as the substrate and the enhanced fluores¬ 
cence due to cleavage of this substrate catalyzed by the 
protease was monitored at 538 nm with an excitation 
wavelength of 355 nm using a fluorescence plate 
reader. The reaction mixture contained 3 pg 3CL pi ° 
protease and 20 pM fluorogenic substrate in 20 mM 
Tris/HCl buffer (pH 7.5) (Grum-Tokars et al. 2008). 
The plates were analyzed at 25°C with continuous 
monitoring of fluorescence for 25 min, with recording 
of relative fluorescence units (RFUs) using a Spectra- 
Max Gemini XPS apparatus (Molecular Devices, 
Sunnyvale, CA, USA) with excitation and fluores¬ 
cence emission wavelengths of 355 and 538 nm, 
respectively. Kinetic parameters of recombinant 
3CL pi ° were obtained using 12.5-100 pM FRET 
peptides in the fluorescent assay with an 18 min 
measurement period. Reaction responses were linear 
within this time. The ( K m ) value was calculated from a 
Lineweaver-Burk using the SigmaPlot program (Sy- 
stat Software, San Diego, CA, USA). 

Inhibition assay 

Quercetin, daidzein, puerarin, epigallocatechin (EGC), 
epigallocatechin gallate (EGCG), and gallocatechin 


gallate (GCG) were purchased from Sigma-Aldrich 
and ampelopsin (AMPLS) was purchased from ZR 
chemicals (Shanghai, China). Primarily inhibitory 
activities of flavonoid compounds were determined by 
measuring the remaining activity of 3CL pi ° at 200 pM 
inhibitors. GCG, EGCG, and EGC were dissolved in 
water; quercetin, puerarin, daidzein, and AMPLS were 
dissolved in dimethylsulfoxide (DMSO) as 10 mM 
stock solutions. The enzyme reaction digest (100 pi) 
was composed of 3 pg enzyme, 20 pM FRET substrate, 
200 pM of each flavonoid compound, and 20 mM Tris/ 
HC1 buffer (pH 7.5). Reactions were run for 18 min at 
25°C with continuous monitoring of fluorescence. 
The inhibition was calculated using following formula 
(1): % inhibition = 100 — remaining activity (%) 
where the remaining activity (%) = [(S — S 0 )f 
0 C - C 0 )] x 100 (1), where C is the fluorescence of 
the control (enzyme, buffer, and substrate) after 18 min 
incubation, C 0 is the fluorescence of the control at time 
zero, S is the fluorescence of the tested samples 
(enzyme, tested sample solution, buffer and substrate) 
after 18 min incubation, and S 0 is the fluorescence of 
the tested samples at time zero. The 50% inhibitory 
concentration (IC 50 ) was defined as the concentration of 
3CL pro inhibitor necessary to reduce 3CL pi ° activity by 
50% relative to a reaction mixture containing 3CL pi ° 
enzyme but no inhibitor. Inhibitor kinetic studies were 
performed for GCG, which was the best inhibitor 
against 3CL pi °. The method was similar to those used in 
a kinetic study of the recombinant enzyme, except for 
the use of multiple concentrations of the inhibitor 
(0-60 pM) and variable concentrations of substrate 
(5-15 pM). The type of inhibition was determined 
using Lineweaver-Burk plots and a Dixon plot (1/v as a 
function of inhibitor concentration, [I]) and kinetic 
parameters (Kj) was calculated using the SigmaPlot 
program. 

Docking of flavonoid compounds with 3CL pi ° 

The three-dimensional structure of 3CL pi ° was 
retrieved from the Protein Data Bank [http://www.pdb. 
org, accession code 2ZU5]. A/-[(benzyloxy)carbonyl]- 
(9-tert-butyl-L-threonyl-A- [(1 R)-4-cyclopropyl-4-oxo-1 - 
{[(3S)-2-oxopyrrolidin-3-yl]methyl} butyl]-L-leucinamide 
(ZU5) was located in the active site of 2ZU5 (Lee et al. 
2009). All water molecules, co-crystal ligand ZU5 were 
removed and the stmcture information containing only the 
amino acid residues of the 3CL pro enzyme was used for 
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docking. Docking files were prepared using AutoDock- 
Tools software (Sanner et al. 1996). For the protein mol¬ 
ecules, polar hydrogen atoms were added and nonpolar 
hydrogen atoms were merged. Kollman charges and sol¬ 
vation parameters were assigned by default. The three- 
dimensional atomic coordinates of GCG were generated 
by the Corina program (Molecular Networks GmbH, 
Erlangen, Germany), Gasteiger charges were added and 
nonpolar hydrogen atoms were merged. The grid box, 

o 

with grid spacing of 0.375 A and dimensions of 
60 x 60 x 60 points along the x, y, and z axes, was 
centered on the macromolecule. AutoDock version 3.0.5 
software using the Lamarckian genetic algorithm (LGA) 
was used for the computational molecular docking sim¬ 
ulation of flexible small molecules to rigid proteins with 
ligand and rigid proteins (Morris et al. 1998). Important 
docking parameters for the Lamarckian genetic algorithm 
were a population size of 250 individuals, maximum of 5 
million energy evaluations, maximum of 27,000 genera¬ 
tions, mutation rate of 0.02, crossover rate of 0.80, and 100 
docking runs (each docking job produced 100 docked 
conformations). The probability of performing a local 
search on an individual in the population was set to 0.06 
and the maximum number of iterations per local search 
was set to 300. The conformation with the lowest docked 
energy was chosen from the most populated cluster and 
was put through to the next stage. The hydrogen bond (H- 
bond) interaction between 3CL pi ° and GCG was identified 
by Ligplot software (Wallace et al. 1995). 

Results and discussion 

Recombinant 3CL pi ° enzyme preparation 

The 918 bp gene encoding 3CL pi ° (amino acids 
3,241-3,546) from human SARS-CoV was cloned 
into the pPICZaA expression vector (pPICZaA- 
3CL pi °) and the pPICZaA-3CL pi ° plasmid was line¬ 
arized by Sad digestion that was further incorporated 
into the AOX1 locus of Pichia pastoris. Seven 
colonies were identified as 3CL pro -positive by PCR 
(Supplementary Fig 2). A single clone was selected 
for expression; the recombinant 3CL pi ° was resolved 
electrophoretically as a band of approx. 42 kDa. To 
confirm that the clone secreted 3CL pi °, western 
blotting was performed using the anti-His antibody; 
3CL pi ° was apparent as the same band upon SDS- 


PAGE and the control did not show this band on 
Western blot (Supplementary Fig 3). 

Large-scale expression of 3CL pi ° was carried out in 
a 10 1 fermenter with 4 1BMMY medium. 3CL pi ° was 
secreted and the amount of 3CL pro protein was 
increased as the induction time increased (Fig. la) as 
confirmed by western blot analysis (Fig. lb). After 
4 days induction, the RFU (as representative of 
enzyme activity) was increased 2.4-times compared 
with 3 days induction. Electrophoretic analysis of the 
3CL pi ° fraction obtained using ammonium sulfate 
gave 80% purity (Supplementary Fig 4). The yields of 
each step in the procedure are summarized in Supple¬ 
mentary Table 1. No appreciable glycosylation was 
observed for 3CL pro even after endoglycosidase H 
treatment. 

To calculate the kinetic parameters of the purified 
3CL pro , enzyme activity was analyzed with fluorescent 
substrate from 12.5 to 100 pM. The Km value was 
15 d= 1 pM (Supplementary Fig 5). The K m value of 
the purified 3CL pi ° was similar to the K m value of 
3CL pi ° expressed in E. coli (K m of 17 ±4 pM) (Kuo 
et al. 2004). 

3CL pi ° inhibition by flavonoid compounds 

Seven flavonoid compounds belonging to four groups 
of flavonoids (flavonol, flavanonol, isoflavone, and 
flavan-3-ol) (Fig. 2) were evaluated for their inhibi¬ 
tory activity against 3CL pi ° expressed from P. pastoris 
GS115. Table 1 shows the inhibitory activity of each 
flavonoid against the 3CL pi ° at 200 pM. Quercetin, 
EGCG and GCG inhibited more than 80% of the 
activity of recombinant 3CL pi °; AMPLS, puerarin, 
and daidzein inhibited more than 30% activity of 
recombinant 3CL pi °; and EGC inhibited only 5% of 
the activity of recombinant 3CL pro . Six compounds 
displayed an IC 50 ranging from 47 to 381 pM 
(Table 1). GCG showed the best inhibition against 
recombinant 3CL pro with IC 50 value of 47 db 0.9 pM, 
and so was used for the further analysis of inhibition 
mode. Both Lineweaver-Burk and Dixon plots were 
used. As shown in Fig. 3a, GCG exhibited competitive 
inhibition toward 3CL pi ° because the Lineweaver- 
Burk plot of 1/v versus 1 /[S] resulted in a family of 
straight lines with the same y-axis intercept. The K x 
value of GCG was determined to be 25 db 1.7 pM 
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Fig. 1 SDS-PAGE and 
Western blot results of 
recombinant 3CL pi ° after 3 
and 4 days induction during 
fermentor culture. Lane M : 
marker; lanes 1 and 2: SDS- 
PAGE after 3 days (lane 1 ) 
and 4 days (lane 2); lanes 3 
and 4\ Western blot 
conducted after 3 days (lane 
4) and 4 days (lane 3 ) 





kDa 

—115,7 
— 79.55 



— 37.13 



_ 29.13 


I 


1 2 M 

from the common v-axis intercept of lines on the 
corresponding Dixon plot (Fig. 3b). 

Molecular docking on 3CL pi ° 

In order to get a better comprehension of the molecular 
recognition process between 3CL pro and antioxidant 
compounds, docking experiments were performed using 
the crystal structure of 3CL pi ° (2ZU5). Autodock 3.0.5 
was used to carry out docking simulations. The free 
binding energy of flavonoid compounds is shown in 
Table 1. Among them, GCG displayed the lowest free 
binding energy (—14 kcal mol -1 ). The binding between 
GCG and active site pocket of 3CL pi ° is shown in 
Fig. 4a. To elucidate the interaction of 3CL pi ° with GCG, 
the potential hydrophobic and H-bond interactions 
between amino acid residues in the active site pockets 
of 3CLpro and GCG were investigated using the Ligplot 
program. Figure 4b depicts the details of the specific 
interactions between GCG and 3CL pi °. Carbon atoms of 
GCG interacted hydrophobically with His41, Cysl45, 
Metl65, Glul66, Aspl87, Argl88, and Glnl89 of 
3CL pro . GCG formed seven hydrogen bonds with 
residues in the catalytic binding pocket of 3CL pi °. The 


— 21.5 


3 4 M 

O atom of the main chain carboxyl group of Glul66 
formed an H-bond with the O atom of the 5-hydroxyl 

o 

group of the A ring with a distance of 2.56 A. The O 
atom of hydroxyl group of aromatic side chain of Tyr54 
formed a H-bond with the O atom of the 4 '-hydroxyl 
group of the ring B at 2.7 A. The O 29 atom of the 5'- 
hydroxyl group of ring B accepted a H-bond from the N 
atom of the imidazol group of His41with a distance 
2.72 A. The O 28 atom of the 3"-hydroxyl group of the 
galloyl group accepted a H-bond with the N atom from 

o 

the imidazol group of His 163 with a distance of 2.96 A. 
The 0“ atom of the 4 -hydroxyl group of the galloyl 
group has two H-bonds: one with the carboxyl group of 
Leu 141 and another one with the carboxyl group of 
Serl44 at 2.90 and 2.49 A, respectively. The N atom of 
the amino group of Glyl43 donated a H-bond with an 
O atom of the 5"-hydroxyl group of the galloyl group 

o 

with a distance of 3.35 A. 

Structural activity relationships of flavonoid 
compounds 

In this study, we compared the inhibition activity of 
AMPLS, EGC, EGCG, and GCG containing the same 
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Fig. 2 Molecular structures 
of the seven flavonoids 
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Table 1 Inhibitory activity of flavonoid compounds against 

o rT pro 


Compound 

Inhibition 3 (%) 

IC 50 (mM) 

Docking score 
(kcal mol -1 ) 

AMPLS 

34 

364 db 8.7 

-9.9 

Quercetin 

82 

73 ± 4 

- 10.2 

Puerarin 

33 

381 ± 12.5 

-11.3 

Daidzein 

34 

351 ± 2.9 

- 8.6 

EGC 

5.4 

ND 

-9.3 

EGCG 

85 

13 ±2 

-11.7 

GCG 

91 

47 ± 0.9 

-14.1 


ND not determined, GCG gallocatechin gallate, EGCG 
epigallocatechin gallate, EGC epigallocatechin, AMPLS 
ampelopsin 

a Inhibition by 200 pM 

B-ring at 200 jiM (Table 1). The decreasing order 
of the inhibitory activity was EGC < AMPLS < 
EGCG < GCG. EGCG and GCG have a galloyl 
moiety at the 3-OH position, which is absent in the 
other catechins used in this study. EGCG and GCG 



EGCG GCG 

displayed stronger 3CL pro inhibitory activity than 
those of EGC and AMPLS. GCG (2S, 3R type), which 
is a C-2 epimeric isomer of EGCG (2R, 3R type), 
showed 1.5-times higher 3CL pi ° inhibitory activity 
than that of EGCG. Molecular docking simulation was 
used to calculate the binding of GCG to the 3CL pi ° 
active site. GCG bound at the substrate-binding pocket 
of 3CL pi ° with numerous hydrophobic and hydrogen 
bond interactions. The galloyl group from GCG was 
important for GCG binding to 3CL pro active site 
pocket because it has four hydrogen bond interactions 
with Leul41, Glyl43, Serl44, and Hisl63 (Fig. 4b). 
The effect of B-ring and hydroxyl group substitution 
on the B-ring for the inhibitory activity was evaluated. 
Daidzein and puerarin, which lack the B-ring, showed 
little inhibition of 3CL pi °. The 3CL pi ° inhibitory 
activity of AMLSP was 4.96-times lower than that of 
quercetin. This confirmed that the addition of an OH 
group at S'-position of the B ring decreased the 3CL pi ° 
inhibitory activity. The effect of the structures of the A 
and C rings on the inhibitory activity was evaluated. 
Since AMPLS, which lacks the 2,3-double bonds in 
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Fig. 3 Lineweaver-Burk plot (a) and Dixon plot (b) analyses 
for the inhibition of 3CL pi ° by GCG. The kinetic constants, K m 
and K h were calculated using linear regression analysis, a GCG 
concentration 0 pM (filled circle ), 20 pM ( open circle ), 30 pM 



[I] &jM) 


(filled inverted triangle), 40 pM ( open triangle ), 50 pM (filled 
square ), 60 pM {open square), b FRET substrate concentrations 
5 pM (filled circle), 7.5 pM {open circle), 10 pM (filled 
inverted triangle), and 15 pM {triangle) 
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Fig. 4 Computational docking and hydrophobic and hydrogen 
bond interactions of GCG with amino acid residues in the active 
site of 3CL pi °. a Comparison of binding modes of GCG {green) 
in the active site pocket of 3CL pi °. b Hydrophobic and H-bond 


interactions between GCG and amino acid residues in the active 
site of 3CL pro . H-bond interactions are represented by green 
dashed lines {Red, oxygen; cornflower blue, nitrogen; black, 
carbon) 
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the C-ring, showed lower inhibitory activity than that 
of quercetin, 2,3-double bonds are probably crucially 
influential to the inhibitory activity. EGC lacking the 
C(4)=0 in the C-ring, C(2)=C(3), containing S'-OH 
group and lacking galloyl moiety showed the lowest 
inhibitory activity compared to that of AMPLS, 
daidzein, puerarin, quercetin, EGCG, or GCG. 

Conclusions 

An important goal of the present study was to 
understand the inhibition by seven flavonoid com¬ 
pounds belonging to four groups of flavonoids (flavo- 
nol, flavanonol, isoflavone, and flavan-3-ol) against 
3CL pro of SARS-CoV. The active extracellular 3CL pro 
was successfully expressed and purified in P. pastoris 
GS115. Among the investigated flavonoids, GCG was 
the best inhibitor against 3CL pi ° by an in vitro assay. 
The structure and inhibition activity relationship 
among seven flavonoid compounds was also investi¬ 
gated. GCG showed numerous hydrophobic and 
H-bonds interaction with amino acid residues in the 
active site pocket of 3CL pro . 
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